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Abstract . In a companion paper, we present observational 
evidence that the stable, growth-phase auroral arc that intensi- 
fies at substorm expansion phase onset often forms on mag- 
netic field lines that map to within ~1 to 2 Re of synchronous. 
The equatorial plasma pressure is 1 to 10 nPa in this region, 
which can give a cross-tail current >0.1 A/m. In this paper, 
we propose that the arc is formed by a perpendicular magneto- 
spheric-current divergence that results from a strong dawn-to- 
dusk directed pressure gradient in the vicinity of magnetic 
midnight We estimate that the current divergence is suffi- 
ciently strong that a > 1 kV field-aligned potential drop is 
required to maintain ionospheric-current continuity. We sug- 
gest that the azimuthal pressure gradient results from proton 
drifts in the vicinity of synchronous orbit that are directed 
nearly parallel to the cross-tail electric field 

Introduction 

A fundamental feature of an auroral substorm is the inten- 
sification and break-up of a pre-existing, stable auroral arc at 
the onset of the expansion phase [e.g., Akasofu, 1964; 
Fukunishi, 1975]. The intensification is accompanied by the 
abrupt decrease in the magnitude of the horizontal magnetic 
field measured on the ground that is traditionally used to iden- 
tify the onset of the substorm expansion phase. 
Understanding the auroral intensification is critical to under- 
standing the physics of magnetospheric substorms. A pre- 
requisite for this is a determination of the processes that lead to 
the formation of the arc that intensifies. However, there has 
been remarkably little attention given to this topic. Here we 
use observed features of the pre-intensification arc and of the 
magnetosphere prior to substorm onset (during the period 
referred to as the "substorm growth phase") to construct a the- 
ory for the formation of the arc. 

The visible arc is extended in longitude and narrow in lati- 
tude (-1 to 10 km). Such narrow arcs typically result from 
precipitating electrons having distributions of the form associ- 
ated with significant (more than a few hundred eV) accelera- 
tion by magnetic field-aligned potential drops, <j)|j [Lyons, 
1992, and references therein]. We presume that the pre-inten- 
sification arc results from such accelerated electrons, so that 
understanding the arc's formation requires understanding the 
formation of a significant <f> ( |. Here we put forward a proposal 
for how this $ | ( is formed over a longitudinally extended 
region. We do not, however, consider how the visible arc 
becomes very narrow. 
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Observational Basis 

In the companion paper [Samson et al., 1992b], we show 
that the arc forms within a 1° to 2° wide latitude region of 
intense atmospheric emissions from proton precipitation. We 
present evidence that the precipitation occurs along magnetic 
field lines that are ~5° equatorward of the open-closed field- 
line boundary and cross the equator at radial distances of ~5 to 
10 R e on the night side. The 1° to 2° latitudinal width for the 
emissions is an upper limit to that of the actual proton precipi- 
tation from the magnetosphere since charge exchange spreads 
emissions from a very narrow beam over nearly 2° in latitude 
[Davidson, 1965]. Mapping this width to the equator implies 
that the equatorial source region for the intense atmospheric 
emissions is <, 1 to 2 Re in radial extent. It has also been 
found [Baumjohann et al., 1981; Robinson and Vondrak, 
1990; Samson et al., 1992a] that the pre-intensification arc lies 
equatorward of the ionospheric reversal in the north-south 
electric field known as the Harang discontinuity. 

Auroral arcs nearly always occur in regions of isotropic (in 
pitch angle) proton precipitation [Lyons et al., 1988], and such 
precipitation is expected to result from violation of the guiding- 
center approximation within the dawn-to-dusk, cross- tail cur- 
rent sheet [Sergeev et al., 1983; Lyons and Speiser, 1982]. 
Sergeev et al.'s calculations show that only modest equatorial 
magnetic-field distortions within the current sheet are neces- 
sary for the precipitation to be isotropic, and that isotropic 
precipitation should occur from the entire radial extent of the 
cross-tail current. Samson et al.'s [1992b] observations show 
that the pre-intensification arc forms in the region of strongest 
precipitating proton energy fluxes. Thus, its equatorial source 
region must be the portion of the cross-tail current sheet that 
has the highest particle pressures during the growth phase. 

The nightside distribution of plasma sheet pressure versus 
radial distance for periods shortly preceding substorm expan- 
sion phase onsets has recently been analyzed by Kistler et al. 
[1992]. Their observations were somewhat off the equator, so 
they presented the total pressure (plasma plus magnetic). They 
argued that the particle pressure should dominate the magnetic 
pressure at the equator, so that the total pressure should be a 
good indicator of the plasma pressure at the equatorial crossing 
of a field line. We adopt this assumption to estimate equatorial 
pressures in the present analysis, though it may lead to some 
errors at the lower radial distances (~5 to 7 R e ) considered here 
due to finite equatorial magnetic pressure. 

Kistler et al [1992] found that the pre-onset pressure 
increases monotonically with decreasing radial distance, 
reaching 1.5 nPa at 10 R e and 4 nPa at 7 R e . They had no 
observations further earthward. Spence et al. [1989] compiled 
observations that show the plasma pressure increasing to ~10 
nPa near 5 Re, though these observations were not sorted by 
substorm phase. Based on these studies, we conclude that the 
highest plasma sheet pressures during a substorm growth 
phase are at equatorial radial distances of -5 to 10 R e . This 
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agrees with the location of the source region for the intense 
proton precipitation inferred by Samson et al. [1992b]. The 
~1 to 10 nPa pressures within this region and their monotonic 
increase with decreasing radial distance presumably result 
from the earthward convection of ions into regions of increas- 
ing magnetic field strength, B. The observation of magnetic 
field distortions near synchronous orbit during the growth 
phase of substorms [McPherron et aL, 1973; Thomas and 
Hedgecock, 1975] is consistent with having high plasma pres- 
sures and with the isotropic precipitation of protons from this 
region via violation of the guiding center approximation. 

To estimate the growth-phase, cross-tail current near syn- 
chronous orbit, we assume pressure balance and a total pres- 
sure of 5 nPa. This gives a magnitude for the x-component of 
the lobe magnetic field of 1 12 nT. (Here x is normal to the 
cross-tail current and parallel to the plane of the tail current 
sheet, and z is normal to the plane of the current sheet). In 
agreement with this, increases in the x-component of B to 
~1 10 nT have been seen during the growth phase of substorms 
at synchronous orbit [McPherron and Manka, 1985, and two 
examples in Figure 7 of Nagai, 1987]. Neglecting the varia- 
tion of B z with x, a lobe field strength of 1 10 nT implies a 
cross- tail current (integrated over z) of magnitude 0. 175 A/m. 
This is consistent with the magnitudes of 0.3 A/m and >0.1 
A/m inferred from magnetic field observations near syn- 
chronous orbit during the substorm growth phase by 
Pulkkinen et al. [1991] and Kaufmann [1987], respectively. 

Theory 

Here we show that a growth-phase plasma pressure of ~ 5 
nPa near synchronous orbit can lead to the formation of the 
pre-intensification arc provided the pressure has a sufficiently 
large azimuthal gradient in the dawn-to-dusk direction. We 
also propose that the required azimuthal gradient exists near 
synchronous orbit and that it results from the predominant 
particle drift in the direction of the cross-tail electric field. 
Here we use "near synchronous orbit" to mean a region within 
about 1 to 2 Re of synchronous orbit. 

The precipitating electrons that form auroras carry an 
upward field- aligned current. However, the maximum current 
density j |, that can be carried by these electrons is given by 
j|l max = en(Kth/27cm e ) 1/2 for <t> 1{ = 0. Here n and K t h are, 
respectively, the density and thermal energy of the magneto- 
spheric electron distribution. For n = 1 cm J and K t h = 300 
eV, j njnax = 5 X 10 -7 A/m 2 . It is necessary that e<t>|, > K t h for 
electrons to carry a larger j ( |. Such a $ ^ increases the electron 
flux within the loss cone for a given magnetospheric value of^n 
and K t h- Auroral arc field lines have j„ > 5 x 10”' A/m , 
which requires that the electrons are accelerated by a <j) (1 > 300 
V [see Lyons, 1992, and references therein]. Thus, under- 
standing the formation of the pre-intensification arc amounts to 
understanding what drives an upward jj, > 5 x 10 _/ a/m . 

To evaluate j || , we make use of the relation [Vasyliunas, 
1970; Birmingham, 1992] between j |( , the plasma pressure P, 
and the flux tube volume V. 

j(, = (Bi/B E )B- VP X VV (1) 

for isotropic pressure. This relation is obtained from the bal- 
ance between j u and the divergence of the magnetospheric cur- 
rent perpendicular to B. Here jj, is the upward current density 


at the ionosphere and Be and Bj are the magnetic field intensi- 
ties at the equator and in the ionosphere; the gradient operator 
is evaluated at the equator, and V = Ids/B is evaluated from 
the equator to the ionosphere. 

Near synchronous orbit, W is directed nearly radially 
outward, and we expect the radial component of VP to be 
directed inward. Thus, it is necessary that VP and/or W 
have an azimuthal component directed from dawn to dusk in 
order to drive an upward field- aligned current. Such an 
azimuthal gradient should also modify magnetospheric electric 
fields and thus particle drifts. This is because j u must be con- 
sistent with current continuity in the ionosphere as well as the 
magnetosphere, and it is necessary to have V«E < 0 in the 
magnetosphere to have an upward directed j n from the iono- 
sphere (Lyons [1992], and references therein). As a result, a 
divergence of the perpendicular magnetospheric current must 
modify the magnetospheric electric field (via small accumula- 
tion of negative charge density) sufficiently to give V*E < 0 
and current continuity in both regions. Iijima et al. [1992] 
suggested that growth-phase, field-aligned currents could 
result from an inward radial component to W. We do not 
consider this possibility here since the increase in plasma pres- 
sure with decreasing radial distance is consistent with the usual 
increase in equatorial B and thus an outwardly directed V V. 

To illustrate our proposal for forming the azimuthal pres- 
sure gradient, we show sketches of equatorial proton trajecto- 
ries and equipotentials in Figure 1. The trajectories are based 
on calculations by Chen et al. [1992] for equatorially mirroring 
20 MeV/Gauss protons (which corresponds approximately to 
the precipitating ion energies in Samson et al. [1992b]), and 
the equipotentials are based on those for a uniform convection 
electric field added to corotation. Both the trajectories and the 
equipotentials have been qualitatively adjusted for expected 
effects of azimuthal pressure gradients. 

Two related features of the ion motion on the night side 
that can give rise to a dawn-to-dusk directed pressure gradient 


■— Source Region for r 71 Harang 

Break-up Arc \Z.A Discontinuity 



Fig. 1. Sketch of proton trajectories and equipotentials in 
equatorial plane during a substorm growth phase. Trajectories 
are based on those for 20 MeV/Gauss protons; equipotentials 
are based on those for a uniform cross-tail electric field plus 
corotation. Both have been qualitatively adjusted for the ex- 
pected effects of azimuthal pressure gradients near the Harang 
discontinuity and the source region for the pre-break-up arc. 



are included in Figure 1. First, as ions convect earthward 
from the tail, magnetic drift and a reduced particle source from 
the tail flanks (relative to that from the far tail) causes a deple- 
tion of ions on the dawn side relative to the dusk side and a 
dawn-to-dusk gradient in particle pressure [Erickson et al., 
1991; Spence and Kivelson, 1992]. This causes a divergence 
of current across the tail that must be balanced by upward 
field-aligned currents. This requires V*E < 0 in the tail, which 
results in an x-component of electric field that reverses in 
direction. Erickson et al. [1991] and Liu and Rostoker [1991] 
proposed that this feature of the ion motion is responsible for 
the Harang discontinuity. The model of Erickson et al. shows 
the electric field reversal to occur at a radial distance of ~10 R e - 

The negatively charged region associated with the Harang 
discontinuity is identified by the shading and associated dis- 
tortion of equip otentials in Figure 1. The x-component of 
electric field in the tail in this region maps to the ionosphere as 
the north-south electric field that reverses direction at the 
Harang discontinuity. The x-component of electric field also 
modifies magnetospheric particle drifts. The earthward- 
directed electric field beyond -10 R e causes a drift towards 
dawn that partially cancels the magnetic drift towards dusk. 
This may allow protons to drift closer to the Earth than they 
would otherwise, which could increase the plasma pressure 
and the cross-tail current near and within 10 Re. 

The second feature of the ion motion that can cause a 
dawn-to-dusk pressure gradient occurs as ions drift earthward 
of 10 Re. Here the magnetic drift speed increases significantly 
relative to the electric drift speed, and ion drifts become pre- 
dominantly azimuthal. The tail ward-directed electric field 
component earthward of 10 R e adds to this azimuthal drift. 
Near midnight, the azimuthal drift causes proton trajectories to 
become approximately parallel to the cross-tail electric field. 
This is illustrated in Figure 1 as occurring between about 5 and 
8 R«. Asa result, protons should gain considerable energy as 
they drift towards dusk, which may lead to a significant 
azimuthal pressure gradient earthward of that considered by 
Erickson et al. [1991] and Liu and Rostoker [1991]. 

We propose that this azimuthal pressure gradient develops 
in the region of large plasma sheet pressure and cross-tail cur- 
rent near synchronous orbit and causes the formation of the 
pre-intensification arc. This proposal is supported by the ob- 
servation [Baumjohann et al., 1981; Robinson and Vondrak, 
1990; Samson et al., 1992a] that the arc lies equator ward of 
the ionospheric mapping of the Harang discontinuity in a re- 
gion of northwest-directed electric fields. The northward 
component is expected from the ionospheric mapping of the 
tail ward-directed component within -10 Re, and the westward 
component is expected from the mapping of the cross-tail 
electric field. Robinson and Vondrak found the westward 
electric field to be 10 to 20 mV/m within the arc. This gives a 
7 to 14 keV energy gain for each hour of local time that a pro- 
ton drifts along the arc’s direction. Kistler et al.’s [1992] ob- 
servations show thermal energies of -20 keV near syn- 
chronous orbit, so that an azimuthal drift of 1 to 3 hr in local 
time should roughly double the energy of a typical proton. 

To evaluate our proposal's feasibility, we use Eq. (1) to 
estimate the intensity of j N and to thus determine whether a 
significant c() i, is required. Kistler et al.'s [1992] observations 
show the equatorial pressure gradient near 7 R e to be dP/dr ~ 
-2 nPa/R^. It would be desirable to evaluate W in a mag- 
netic field model self-consistent with an assumed pressure 
distribution. However, here we obtain an estimate for dV/dr 
by simply taking its dipole value 1.83L /B 0 ur/Re, where B 0 
is the equatorial surface-field strength. To estimate the az- 


imuthal gradient of P, we assume that the pressure increases 
by 2nPa over 2 hr in local time. We also assume that the mag- 
netic field stretches as P increases so that the equatorial mag- 
netic field strength decreases as P increases. This gives an 
azimuthal gradient in V , which we take to be a doubling of 
the dipole value over 2 hr in local time. (We cannot find much 
literature on nightside azimuthal pressure or magnetic field 
gradients, but Fairfield [1986] found equatorial magnetic field 
magnitudes on the dusk side of the tail to be about 1/3 of those 
near midnight at -20 < x < —10 R<0 For L = 7, these values 
give the estimate that j ^ = 3 x 10° A/m 2 , the contributions 
from the assumed azimuthal gradients of P and V being ap- 
proximately equal. This is 6 times greater than the maximum 
current density estimated above for = 0, and it thus re- 
quires a <{> || of a few kV. This is required to form the pre- 
intensification arc. 

The distortion of equipotentials in the vicinity of the arc 
source region is schematically illustrated in Figure 1 . In the 
figure, we show the source region for an actual arc to be nar- 
rower than the 1 to 2 Re region of distorted potentials associ- 
ated with the upward field-aligned current; however, we offer 
no explanation at this time for this nanowness. 

Conclusions 

We have discussed observational evidence that the auroral 
arc that intensifies at substorm onset is formed on magnetic 
field lines that map to within -1 to 2 R e of synchronous orbit, 
where plasma pressures are -1 to 10 nPa during the substorm 
growth phase. These large pressures cause an enhanced 
cross-tail current, which may have a magnitude >0.1 A/m. 
This current is associated with a distortion of the magnetic 
field that should be sufficient to isotropize proton pitch-angle 
distributions via violation of the guiding center motion. This 
will cause the proton precipitation observed in the vicinity of 
the pre-intensification arc. 

We propose that the plasma pressure in this region has a 
strong dawn-to-dusk directed gradient, which drives a field- 
aligned current upward from the ionosphere. We estimate that 
the magnitude of jj ( is sufficiently large that a field-aligned 
potential drop 4> ( | > 1 kV is necessary for precipitating elec- 
trons to carry the current. This potential drop energizes the 
precipitating electrons that form the arc. As illustrated in 
Figure 1, we propose that the pressure gradient results from 
the azimuthal proton drift near synchronous orbit, which is in 
the direction of the cross-tail electric field. 

We expect the pre-intensification arc to form in the region 
of northwest-directed (in the northern hemisphere) electric 
field equatorward of the Harang discontinuity, as has been 
observed. Also, to maintain current continuity in the iono- 
sphere, we expect the northward component of the electric 
field to be reduced in the vicinity of the arc, as observed by 
Robinson and Vondrak [1990] and illustrated in their Figure 
14a (also see Samson et al., [1992a]). The related electric 
potential distribution in the equatorial plane is illustrated in 
Figure 1. It is reasonable to expect that the substorm expan- 
sion phase is initiated on the magnetic field lines of the arc that 
intensifies. Under this assumption, our proposal is consistent 
with recent observational evidence [e.g., Lopez et al., 1990, 
Lui et al., 1992] that a substorm expansion phase is initiated 
near synchronous orbit via cross-tail current disruption. 

Finally, since we expect the protons in the vicinity of syn- 
chronous orbit to be precipitated with approximately isotropic 
distributions, it may be possible to relate proton precipitation 
measurements to the equatorial plasma pressure during the 
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growth phase. For an isotropic equatorial particle distribution, 
pressure P = (47t/3)mjfv 4 dv, and energy flux per steradian £ = 
(m/2)Jfv 5 dv, where m is mass and f is the particle distribution 
function. For a Maxwellian distribution of thermal speed v^, 
we obtain £ = (v t h/rc 3 ^)P. The total precipitating energy flux, 
£p = 27cj£cosasinada = 7t£, where a is pitch angle, so that 
for P = 1 to 10 nPa and a proton thermal energy of 20 keV, we 
obtain £p = 10” 3 to 10” joules/m 2 -s. The optical observa- 
tions of Fukunishi [1975] and Samson et al. [1992b] imply a 
proton energy deposition into the atmosphere of 10 
joules/m 2 -s (based on 70 R of Hp emission per 10 -3 
joules/m 2 -s [Strickland, private communication, 1992]). This 
represents a minimum estimate for £ p because of the ~2° 
spread in the emissions that results from charge exchange. 
Proton measurements from polar-orbiting Stellites could be 
used to search for the large values of £p we expect. However, 
such measurement would need to extend to significantly higher 
energies than are typical. For example, a 20 keV Maxwellian 
would have 70% of £p at energies above the 30 keV upper 
limit on DMSP satellites, and the fraction would be even 
higher if a high-energy tail to the distribution were included. 

Acknowledgments . This work was supported in part by 
NSF grant ATM- 17637 and NASA grant NAG 5-1541 (Space 
Physics Theory Program) to The Aerospace Corporation. 
Research by J. C. Samson was supported by the Natural 
Sciences and Engineering Council of Canada. 


References 

Akasofu, S.-I., The development of the auroral substorm, 
Planet . Space Sci ., 12, 273, 1964. 

Baumjohann, W. J. Pellinen, H. J. Opgenoorth, and E. 
Nielsen, Current systems associated with local auroral 
break-ups. Planet . Space Sci., 29, 431, 1981. 

Birmingham, T. J., Birkeland currents in an anisotropic, mag- 
netostatic plasma, /. Geophys. Res., 97, 3907, 1992. 

Chen, M. W., M. Schulz, L. R. Lyons, and D. J. Gomey, 
Stormtime transport of ring-current ions, J. Geophys . 
Res., 1992 (submitted). 

Davidson, G. T., Expected spatial distribution of low-energy 
protons precipitated in the auroral zones, J . Geophys. 
Res., 70, 1061, 1965. 

Erickson, G. M., R. W. Spiro, and R. A. Wolf, Physics of 
Harang discontinuity, J. Geophys. Res., 96, 1633, 1991. 

Fairfield, D. H., The magnetic field of the equatorial magneto- 
tail from 10 to 40 R^, J. Geophys. Res., 91, 4238, 1986. 

Fukunishi, H., Dynamic relation between proton and electron 
auroral substorms, J. Geophys. Res., 80, 553, 1975. 

Iijima, T. et al., Substorm currents in the equatorial magneto- 
sphere, J. Geophys . Res., 1992 (submitted). 

Kaufmann, R. L., Substorm currents: growth phase and 
onset, J. Geophys. Res., 92, 7471, 1987. 

Kistler, L. M. et al.. Pressure changes in plasma sheet during 
substorm injections,/. Geophys. Res., 97, 2973, 1992. 

Liu, W. W., and G. Rostoker, Effects of dawn-dusk pressure 
asymmetry on convection in the central plasma sheet, /. 
Geophys. Res., 96, 11501, 1991. 


Lopez, R. E. et al., Multipoint observations of a small sub- 
storm, /. Geophys. Res., 95, 18,897, 1990. 

Lui, A. T. Y. et al.. Current disruptions in the near-Earth neu- 
tral sheet region, /. Geophys. Res., 97, 1461, 1992. 

Lyons, L. R., Formation of auroral arcs via magnetosphere- 
ionosphere coupling. Rev. Geophys., 30, 93, 1992. 

Lyons, L. R. and T. W. Speiser, Evidence for current-sheet 
acceleration in the tail,/. Geophys. Res. 87, 2276, 1982. 

Lyons, L. R., J. F. Fennell and A. L. Vampola, A general 
association between discrete auroras and ion precipitation 
from the tail, /. Geophys . Res., 93, 12,932, 1988. 

McPherron, R. L. and R. H. Manka, Dynamics of the 1054 
UT March 22, 1979, substorm event: CDAW 6, /. 
Geophys. Res., 90, 1175, 1985. 

McPherron, R. L., M. P. Aubry, C. T. Russell, and P. J. 
Coleman, Satellite studies of magnetospheric substorms 
on August 15, 1968, /. Geophys. Res., 80, 3068, 1973. 

Pulkkinen, T. T. et al.. Modeling the substorm growth phase 
using the Tsyganenko model and multi-spacecraft obser- 
vations: CDAW-9, Geophys. Res. Lett., 18, 1963, 1991. 

Robinson, R. M., and R. R. Vondrak, Electrodynamical 
properties of auroral surges, /. Geophys. Res., 95, 7819, 
1990. 

Samson, J. C. et al., Substorm intensifications and field line 
resonances in the nightside magnetosphere, /. Geophys. 
Res., 8495, 1992a. 

Samson, J. C. et al., Proton aurora and substorm intensifica- 
tions, Geophys. Res. Lett., 1992b (this issue). 

Sergeev, V. A. et al., Pitch-angle scattering of energetic pro- 
tons in the magnetotail current sheet as the dominant 
source of their isotropic precipitation into the nightside 
ionosphere. Planet. Space Sci., 31, 1147, 1983. 

Spence, H. E., M. G. Kivelson, R. J. Walker, and D. J. 
McComas, Magnetospheric plasma pressures in the mid- 
night meridian: observations from 2.5 to 35 R , /. 
Geophys. Res., 94, 5264, 1989. e 

Spence, H. E. and M. G. Kivelson, Contribudons of the low- 
latitude boundary layer to the finite width magnetotail con- 
vection model, J. Geophys. Res., 1992 (in press). 

Thomas, B. T. and P. Hedgecock, Substorm effects on the 
neutral sheet inside 10 Earth radii, in The 
Magnetospheres of Earth and Jupiter, ed. by V. 
Formisano, pp. 55-70, D Reidel Publishing Co., 
Dordrecht-Holland, 1975. 

Vasyliunas, V. M., Mathematical models of magnetospheric 
convection and its coupling to the ionosphere, in Panicles 
and Fields in the Magnetosphere, ed. by B. McCormac, 
p. 60, Hingham, Mass., 1970. 


L. R. Lyons, Space and Environment Technology Center, 
M2-260, The Aerospace Corporation, P. O. Box 92957, Los 
Angeles, California, USA 90009-2957. 

J. C. Samson, Canadian Network for Space Research, 
Department of Physics, University of Alberta, Edmonton, 
Alberta, Canada, T6G 2J1. 

(Received July 27, 1992; 
revised September 25, 1992; 
accepted October 2, 1992. 


4 


